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Abstract. - The statistical properties of the conductance of one dimensional disordered 
systems are studied at finite bias voltage V and temperature T, in an independent-electron 
picture. We calculate the complete distribution of the conductance P{G) in different regimes 
of V, T within a statistical model of resonant tunneling transmission. We find that P{G) 
changes from the well-known log-normal distribution at T = in the linear response regime to 
a Gaussian distribution at large V, T. The dependence on T and V of average quantities such as 
(G), (InG) is analyzed as well. Our analytical results are confirmed by numerical simulations. 
We also discuss the limits of validity of the model and conclude that the effects of finite T, V 
presented here should be observable. 



Quantum electronic transport in one-dimensional disordered systems has been widely stud- 
ied ever since the concept of Anderson localization has been introduced in the sixties [1]. In 
particular it is known that within the model of noninteracting electrons all states are localized 
in Id, with average extension ^, the localization length. For systems of finite length L, the 
typical conductance decreases exponentially with increasing length. However, there are large 
sample to sample fluctuations and a statistical analysis of the conductance G is required. For 
a Id wire of length L ^ ^, the full distribution of linear conductance at T = has been 
determined [1-5]: it is log- normal with a large variance of the order of the mean. 

Finite temperature T and a finite bias voltage V have a dramatic effect on the conductance 
distribution, as we show in this letter. The conductance is given by an energy integral of the 
transmission coefficient g{E). A window of energies within which electron transport is possible 
is determined by the value of T and V. Since g{E) is a strongly fluctuating function of E, the 
result depends sensitively on the width of the energy window. 

In the linear conductance regime (zero bias voltage), Azbel and Vincenzo [6] calculated 
the T dependence of the averages of the conductance and its logarithm, at low temperatures. 
They employed a model of g{E) as a sum of resonances with vanishing width (J-functions) 
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caused by resonance tunneling through locahzed states situated at random positions, and 
decaying exponentially over the distance ^. They did not consider the full distribution of 
conductances. Moreover, they did not take into account the fluctuations of the localization 
length. Their model therefore does not recover the log-normal distribution of G{T) at T = 0. 
More recently, it has been shown numerically that the average resistance decreases strongly 
with increasing temperature [7] . 

The statistical properties of the nonlinear conductance G{T, V) = J{T, V)/V, in particular 
its distribution P{G) have not been considered before (here J(T, V) is the charge current), to 
our knowledge. 

One may ask whether interaction effects will preempt the behavior found here within the 
independent particle model. The phase coherence required for resonant tunneling is limited to 
a finite length by phase relaxation at finite tcmpcraturc/voltagc. Wc will argue below that 
for sufficiently short wires L L^, even at the values of T, V of interest and phase coherence 
is preserved. Furthermore, at finite T, V transport will occur by way of inelastic processes, 
e.g. thermally activated hopping [8-10]. As shown below, these processes will dominate the 
transport in a regime of T, V beyond the one we will consider here. Our estimates show that 
the model we consider should be applicable to realistic systems. 

In this letter wc study the effect of bias voltage and temperature on the statistical prop- 
erties of transport in a one-dimensional disordered system, within the model of independent 
electrons. We consider transmission through the wire by way of resonant tunneling, and 
employ a statistical model of individual resonances with Brcit-Wigncr line shapes [11]. Wc 
calculate the full distribution of the conductance in various regimes of T, V, as well as average 
quantities such as (G) and (InG). 

We show that the distribution of the conductance develops from a log-normal distribution 
in the limit of zero voltage and temperature to a Gaussian distribution in the high V and T 
regime. We also find that the conductance average {G) is independent of V,T, while (— InG) 
decreases with V, T from its zero temperature and voltage value up to — ln(G), at high values 
of V, T. 

Our analytical results are supported by numerical simulations of a disordered one- dimensional 
multibarrier system. We calculate numerically the total scattering matrix of the system by 
combining the scattering matrices associated to each barrier. We consider all the barriers to 
have the same fixed height, but random widths and separations [12]. 

The model. - The nonlinear conductance G{T,V) in units of Go = 2e^/h is given by 
[13,14] 

where f{E) = (exp[(£ — /i)/fcT] + 1)~^ is the Fermi function; k, the Boltzmann constant; and 
/U, the chemical potential. g{E) is the transmission probability, equal to the zero temperature 
linear conductance for E = Ep = /x(T = 0). In a wire of length L :§> ^, where ^ is the average 
localization length, transmission through the wire takes place via localized states with random 
energies E,y, located at random positions in the wire. The wave function of any of these 
localized states decays exponentially over a random length 

The transmission probability associated with each localized state is given by a Breit-Wigner 
formula, with width T,, = vi^^ + Ti''\ T^f'''^ = Aexp[-(i ± 2z^)/^^]. vi^'''^ are the partial 
widths induced by coupling to the left/right lead, and A is the mean- level spacing of resonant 
states. We will assume <c: A in the following, corresponding to L/^ » 1. Within this 
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model, the conductance g{E) is given by 

On resonance {E = E^), the terms in the sum eq. (0), ijy, depend on the location of the state, 

= [2cosh(2z,/^,)]-2 (3) 

and are maximum for = (center of wire). 

We assume Ei, and z^, to be uniformly distributed in the intervals \z\ < L/2 and [v — 
< E,y < {i^ + 5)A, respectively. The inverse localization lengths (more conveniently 
Xi, = 2L/^^) are assumed Gaussian distributed 

p{x^) = Cexp [-{x^ - {x^)f/4{x^)], (4) 

where = 2L/S, and C is a normalization constant. 

The above model of a disordered Id system leads to the well-known log-normal distribution 
for the linear conductance at T = and L 3> ^ [15]: 

P(ln.g) =C'exp[-(l/2)(4i/e)-i(ln(l/5)-2L/e)2] ; (5) 

C" is a normalization constant. In this case the resonance closest to the Fermi level, Ep, 
dominates and one may safely neglect additional resonances. 
At finite T, V, the nonlinear conductance is given by 



G = -27rfcTRe 



EE 



t r 



,T^o {E. - /i + + zc^„)2 - (^) _ 



(6) 



where we have performed the integrals in eq. by using the expansion 
djjdE = 2kTRc{J2n=aiE " M " «^^n)"^} with w„ = (2n + l)nkT . 

Gaussian distribution at high temperatures/voltages. - Let us first consider the regime 
of not too low temperatures defined by ^ kT. Neglecting in the denominator of eq. 
© yields 

f F 



eV 



_^ eV eV 

tanh — tanh — — 

2kT 2kT 



(7) 



where E^, = E^ — fi. 

Now, for kT or eV ^ A, many terms in the sum over with < i?max, will contribute, 
where -Bmax = coth(ey/2fcT) ~ max(ey, 2fcT). Consider E'max = eV, for example. In this 
case, G(T, V) is given by 

G(T, « I A ^ gME.u.. - E,), (8) 

here gi, — t^F^/A and &{x) is the step function. Since the quantities g^, are statistically 
independent, we may apply the central limit theorem to find a Gaussian distribution for the 
conductance : 

P(G) = Cexp [-(G -(G)) W], (9) 
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Fig. 1 - Distribution of the linear conductance at high temperatures. The solid line correspond 
to a Gaussian distribution with average and variance taken from the numerical simulation. Inset: 
variance vs 1/kT, at high temperatures. A straight line is well fitted to the numerical data, verifying 
the predicted 1/kT behavior of the variance. 



where (G) = (g^) oc cxp (— is essentially the zero temperature and voltage result, while 
the variance = {{gi, — {gy)))'^ decreases as a function of the inverse voltage. C is a 
normalization constant. A similar analysis may be done for the linear conductance regime. In 
this case i?max = kT and the average of the Gaussian distribution is given by {g^) (independent 
of T), while decreases with the temperature as T"'^. In fig. ^we show an example of the 
central limit theorem at work: P{G) from our simulation of a Id system follows a Gaussian 
distribution and its variance (inset) dependent on the temperature as as predicted 

above [16]. 

We have shown that the average conductance at high temperature/voltage is given by 
its corresponding zero temperature result. In fact, our model gives a very smooth tempera- 
ture/voltage dependence for {G{T,V)) at all regimes of T,V. From eq. ©, averaging over 
disorder we find for i ^ ^ 

/nf^ m\ ~ 2^^/ \ 1 cosh wjv^ cosh wi_ 

G(T, V)) « —r^{g)Q In , (10) 

eV coshwAr_ coshwi^ 



where {g)o = •\/^/27rL exp(— L/2^) is the average linear conductance at zero temperature, 
WN+,_ = iA/2kT){N ± eV/2A), and wi^ _ = (A/2fcT)(l ± eV/2A), with N the total num- 
ber of resonances that contributes to the sum in cq. © {N ^ E'^ax/A, at high tempera- 
tures/voltages), eq. (fTn|l depends weakly on T, V: for example, we can see that in both limit 
cases eV/kT » 1 or eV/kT < 1, (G(T, V)) « (5)0. In contrast, (In G(T, V)) depends strongly 
on T,V. At high temperatures/voltages, we have seen that the width of the Gaussian dis- 
tribution P{G) decreases with the inverse temperature/ voltage, keeping constant its average, 
i.e. P{G) becomes a narrow function around (G). As a consequence, (InG) approaches to 
ln(G) at large values of T/V. As a verification of the above results, in fig. |2]we show the 
results from the numerical simulation: the value (G) at different finite temperatures lies on 
the corresponding zero temperature result (solid line). In the inset, (— InG) = 2L/^ at T = 0, 
while (InG) approaches to ln(G) (solid fine), at high T. 
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Fig. 2 - (G(r)) from the numerical simulation as a function of the length of the system (in units 
of ^) at different temperatures. The solid line is the zero temperature result (17)0 (below eq. 1101 . 
Similarly, the inset shows (— InG) at several temperatures, here the solid line correspond to ln{g)o 
(the dotted lines are to guide the eye). 



Intermediate temperatures/voltages. - When T^, <^ E-^aayi £ A, only the resonance closest 
to the Fermi level, say v — will contribute appreciably to G, and 



tanh — " „. „ ^ tanh — ^ 



eV/2 



2kT 



2kT 



The distribution of conductances is then given by 



f-oo pL/2 /•A/2 

P(G) = C / dxp{x) / dz dEaS (G - G(T, V)) , 

Jo J-L/2 J-A/2 



(11) 



(12) 



with p{x) given by eq. G is a normalization constant. Performing the integral over Eq 
with the help of the 5-function we find 



F(G) = G' / dxpix) / ' ds ^° 

•^^1 -^0 A/df § - cotht;)2 - eschew 



(13) 



where go = 2e~^^^csch{xs), v = eV/2kT, Sc = arccosh(6exp(— a;/2))/x, xi = ln(26 — 1), 
X2 = 2 Info, and b = {-K/S./eVG) sinhu/ (1 + coshu). 

In fig. 131 (insets) we compare P(ln G) as calculated from eq. (P(ln G) = GP{G)) to the 
numerical simulation data, in the linear and nonlinear conductance regime. The agreement 
is seen to be good. Compared to the T and V zero result, the peak of P(lnG) is shifted to 
larger values of G (smaller values of InG), and is clearly narrower. 

Typical conductance at low temperatures/voltages. - When i?max ^ Pi/ <IC A, there is 
no significant effect of finite temperatures/ voltages at the level of the bulk distribution of the 
conductance: P(lnG) is seen to be well described by the zero temperature linear conductance 
distribution, as it is shown in Fig^ To discuss the correction to this result we now calculate 
(InG). At small values of V and T wc find, from eq. ((HJ, 



9 T 

-lnG(F,r))« — 



6A^ 



[{eVf + {2nkTf 



(14) 
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Fig. 3 - P(lnG) from numerical simulations (histograms) at 2L/^ — 7. Left and right insets show 
an example of the distribution in the linear and nonlinear conductance regime, respectively. A good 
agreement with the analytical result eq. (I13t is seen (solid line). Main frame: The effect of finite 
T, V on P(lnG) for a case where kT < A < eV (see conclusions). The values of T, V used here are 
beyond the validity of eq. (I13II . in dashed line. As a reference, we also show P(lnG') at zero T,V 
(histograms in dotted line). 



The first term in eq. (|14ll corresponds to the zero temperature result for the linear conductance 
regime, as one may expect. In figure^] (insets) we compare our expression (|14(l to the numerical 
simulations. 

Observability of effects and conclusion. - We now estimate the effect of inelastic processes 
at temperatures kT ^ A. At temperatures kT <^ {Tivp/aO: where a < 1 is a dimensionless 
interaction constant, the phase relaxation length ~ ^ ^ provided 




Fig. 4 - P(G) in the linear regime at kT = 10~^ and 2L/^ ~ 7. The solid line correspond to the 
zero temperature result. Insets: (InG) from the numerical simulation (dots) at small temperatures 
iy = 0) and small voltages (T = 0) . The dashed line is a fit of a function as given by eq. The 
fit parameters (2L/^, A) are consistent with our numerical calculations [16]. 
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kT < ^1 where Ep/l^ is very large; thus for kT-^A phase coherence is seen to 

be preserved. The contribution of variable range hopping has been estimated in [9, 10]. It is 
negligible compared to resonance tunneling at temperatures kT < A(^/L). This means that 
in the regime kT > A thermally activated hopping processes dominate, and the narrowing of 
P(G) in the linear response regime will not be accessible. At kT ^ A and eV ^ A, however, 
thermal activation will not be relevant. Also, the phase coherence should be maintained as 
long as eV <C (■^)"A, which would allow to access the Gaussian regime of P{G). In fig. |3| 
(main frame) we show P(ln G) for values of T, V , as discussed above. A narrower distribution 
is clearly seen, compared to the linear conductance regime. On the other hand, some of the 
effect of inelastic scattering may be included in the present model by adding an inelastic 
component to the Breit-Wigner resonance width [9]. 

To conclude, we expect that the statistical effect of averaging over the resonance structure 
of the transmission coefhcient at large V, T should lead to large observable changes in the 
conductance distributions of quantum wires. Experimental observation of this phenomenon 
should be of considerable interest [18]. 

This work has been supported in part by the DFG Center for Functional Nanostructures 
at the University of Karlsruhe. We thank K. A. Muttalib for helpful discussions. 
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